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I. 


INTRODUCTION 


The  newsboy  problem  was  described  by  Morse  and  Kimball 
early  in  1950  [1].  This  problem  has  been  presented  in  the 
literature  under  a  variety  of  names,  including  newsvendor 
problem,  Christmas  tree  problem,  and  high-fashion  apparel 
problem.  This  kind  of  problem  may  be  described  in  the 
following  way:  The  newsboy  has  to  decide  how  many  papers  to 
purchase  for  resale.  If  he  buys  too  many  papers  he  will 
incur  loss  due  to  left-over  papers  (little  salvage  cost) . 

If,  on  the  other  hand,  he  buys  too  few,  he  should  also  incur 
an  opportunity  cost  due  to  lost  sales.  To  determine  the 
optimal  order  quantity,  it  is  necessary  to  balance  two 
opposing  costs.  Solutions  to  the  newsboy  problem  under  both 
conditions  of  risk  and  uncertainty  are  well  known.  The 
determination  or  estimation  of  the  demand  distribution  is 
critical  in  dealing  with  this  problem.  Optimal  solutions  to 
the  newsboy  problem  which  minimize  expected  cost  require 
information  about  a  particular  quantile  of  the  demand  dis¬ 
tributions.  If  the  demand  follows  some  known  distribution 
the  problem  can  be  solved  by  decision  making  under  risk,  as 
is  also  the  case  when  the  key  quantile  is  known.  With  only 
information  concerning  the  range  of  demand,  optimal  solutions 
can  be  found  by  some  procedure  such  as  the  minimax  rule  under 
uncertainty. 

The  purpose  of  this  study  may  be  summarized  as  follows. 

We  are  concerned  with  decision  policies  during  a  period  when 
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demand  data  is  being  acquired  at  the  rate  of  one  observation 
per  decision  period,  starting  with  no  data  except  an  estimate 
of  the  maximum  demand.  To  find  the  optimal  expected  value 
decision  rule,  five  candidate  quantile  estimators  will  be 
proposed  on  the  basis  of  order  statistics.  They  will  be 
evaluated  by  comparing  their  performance  with  each  other  and 
with  the  well-known  minimax  rule  under  uncertainty.  An  esti¬ 
mator  which  has  the  smallest  relative  cost  deviation  to  the 
risk  cost  solution  (ideal  case)  will  be  proposed  as  optimal. 
We  shall  also  examine  the  switching  time  to  change  from  the 
uncertainty  rule  to  applying  the  proposed  estimator. 

In  Chapter  II,  we  shall  review  the  newsboy  problem  solu¬ 
tions  under  conditions  of  both  risk  and  uncertainty,  and  then 
suggest  alternative  approaches  for  demand  estimation,  and 
for  the  switching  period  from  uncertainty  to  risk.  For  esti¬ 
mation,  five  quantile  estimators  will  be  proposed  in  Chapter 
III.  In  Chapter  IV,  candidate  quantile  estimators  will  be 
evaluated  by  simulation.  The  performance  of  the  estimators, 
and  candidate  switching  periods  will  be  discussed  in  Chapter 
V.  Conclusions  and  recommendations,  for  further  work,  will  be 
given  in  the  final  chapter. 
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II.  STATEMENT  OF  THE  PROBLEM 


In  this  section  we  shall  briefly  describe  the  well-known 
newsboy  problem,  illustrating  its  optimal  solutions  under 
both  risk  and  uncertainty  with  respect  to  demand.  Then  we 
shall  describe  the  quantile  estimation  problem  we  wish  to 
solve . 

A.  THE  NEWSBOY  PROBLEM 

The  problem  may  be  generalized  in  the  following  way. 

The  decision  maker  selects  a  number  (the  newsboy's  order 
quantity)  but  the  actual  demand  (the  number  of  papers  his 
customers  will  want  to  purchase)  may  have  a  randomly  deter¬ 
mined  value.  If  the  decision  maker's  number  is  greater  than 
this  he  pays  a  cost  proportional  to  the  difference,  or  he 
pays  a  cost  proportional  to  the  amount  by  which  his  number  is 
smaller  than  this.  His  cost  is  zero  only  if  his  number  is 
the  same  as  the  actual  demand. 

Let  demand  D  be  a  continuous  random  variable  and  let 
quantity  S  be  the  number  the  decision  maker  selects  to  have 
on  hand.  The  newsboy  cost  equation  may  be  formulated  as  a 
two-piece  continuous  linear  function  in  the  following  manner. 


C(S) 


C  (S  -  D)  if  0  <_  D  _<  S 
s 

< 


„  CQ  (D  -  S)  if  D  >  S, 


(1) 
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where  C  is  unit  cost  of  surplus  when  S  exceeds  D,  and  C 
s  o 

is  unit  cost  of  outage  or  shortage  when  D  exceeds  S.  The 
cost  function  is  shown  in  Figure  1. 

Cost  C(S) 


Figure  1.  The  Newsboy  Cost  Function 

B.  THE  NEWSBOY  PROBLEM  SOLUTIONS  UNDER  RISK 

For  our  case  where  demand  is  a  continuous  random  varia¬ 
ble  D  with  probability  density  function  f (D) ,  and  S  is  a 
decision  variable  whose  value  the  decision  maker  is  trying 
to  find  over  D,  the  expected  cost  may  be  represented  as: 

s  00 

E  [C  (S)  ]  =  C  j  (S-D)f(D)dD  +  C  f  (D-S)f(D)dD. 

s  J  O  ' 

0  3 

(2) 

We  wish  to  find  quantity  S,  which  minimizes  the  expected 
cost.  The  optimum  (minimum  in  this  case)  value  of  S  may 
be  found  by  differentiating  the  expected  cost  function  (2) 
with  respect  to  S,  setting  the  derivative  to  zero,  and  solving 
for  the  optimum  value  S*: 
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dE [C(S) ] 

dS 


oo 


set 


s 

=  C  /  f  (D)  dD  -  C  f  f  (D)  dD  =  0.  (3) 

S  0  °  s 


From  this  we  obtain  the  optimal  value  S*  using  the  cumula¬ 
tive  distribution  function 


X 

F(X)  =  /  f  (D)  dD, 

0 


and  the  result  is  the  well-known  result, 

F(s*>  -  rrr  •  (4) 

s  o 


Differentiating  Equation  (3)  with  respect  to  S  to  get  the 
second  derivative  at  S*, 


cTE[C(S)  ] 


d'S‘ 


S=S’ 


(C+C  )  f(S*)  . 

s  o 


Hence,  S*  is  the  minimum  expected  cost  solution  since 
(C  +C  )f(S*)  is  positive. 

b  U 

To  summarize,  the  optimal  solution  S*  to  minimize  the 
expected  cost  is  the  [CQ/ (C  + CQ) ] th  quantile  of  the  demand 
distribution,  as  shown  in  Figure  2,  and  the  crucial  infor¬ 
mation  we  need  about  the  distribution  is  that  quantile. 

C.  THE  NEWSBOY  PROBLEM  SOLUTIONS  UNDER  UNCERTAINTY 

The  decision  maker  may  be  confronted  with  a  newsboy 
problem  without  knowing  or  being  able  to  estimate  the 
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F(D) 


Figure  2.  Minimum  Expected  Cost  Solution 

probability  distribution  of  demand.  The  problem  of  deciding 

on  S  in  this  situation  is  called  a  decision  under  uncertainty. 

The  decision  maker  may  be  able  to  estimate  an  upper  bound 

for  the  demand  value.  In  this  case  we  say  that  all  we  know 

about  demand  is  that  demand  D  is  such  that  0  <  D  <  D 
—  _  —  —  max 

Three  well  known  approaches  to  decision  making  under 

—  ■"  ~*\ 

uncertainty  are  Laplace  solution  under  assumption  of  a  uni¬ 
form  distribution  of  demand,  minimax  cost  solutions  and 

2, 

minimax  regret  solutions  [2] .  For  the  newsboy  problem, 
these  all  lead  to  the  same  rule  to  find  the  optimal  solution. 
We  shall  discuss  only  the  minimax  cost  approach.  This 
approach  is  to  choose  S  so  that  the  worst  possible  cost  will 
be  minimized,  that  is,  tO'  minimize  maximum  cost.  The  maxi¬ 
mum  cost  occurs  at  D  =  0  or  at  D  =  D  ,  depending  on  the 
cost  coefficients.  This  can  be  seen  from  Figure  1.  If 
maximum  cost  occurs  at  D  =  0 ,  we  would  want  to  reduce  S 
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while,  if  maximum  cost  occurs  at  D  =  D  ,  we  would  want 

max 

to  increase  S  to  reduce  cost.  We  can  minimize  the  maximum 

cost  when  S  is  chosen  so  that  the  cost  at  D  =  0  is  equal  to 

the  cost  at  D  =  D  .  Equating  these,  we  obtain: 

max  ^  ^ 


C 

s 


S* 


C  (D 
o  max 


S*) 


or 


S* 


c 

( _ — 

+  c 

s 


-) 

o 


D 

max 


/ 


(5) 


So  far,  we  have  introduced  the  minimum  expected  cost 

solution  under  both  risk  and  uncertainty  conditions.  The 

minimax  approach  is  tied  to  our  initial  estimate  of  D  , 

max 

and  thus  the  optimal  value  S*  from  the  minimax  rule  will 

be  in  error  to  the  extent  that  our  estimation  of  D  is 

max 

in  error.  This  error  will  continue  over  many  decision 

periods  unless  our  estimate  of  D  is  revised. 

max 

D.  QUANTILE  ESTIMATION  AND  A  DECISION  RULE  SWITCHING 
POINT  FOR  THE  NEWSBOY  PROBLEM 

We  are  interested  in  finding  a  quantile  estimator  which 
gives  better  results  than  the  minimax  approach  under  condi¬ 
tions  of  limited  demand  data,  and  also  in  finding  the 
switching  period  which  is  the  time  at  which  we  switch  from 
the  minimax  rule  to  the  expected  cost  rule  using  the  quantile 
estimator.  The  purposes  of  this  study  are  to  find  the  "best" 
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quantile  estimator  among  five  candidate  estimators,  and  to 
obtain  information  about  the  switching  period. 

We  recall  that  we  wish  to  find  the  quantile  of  the 
cumulative  distribution  function,  F(D).  The  optimal  S*  is 
found  by  an  inverse  mapping  of  the  [CQ/ (Cg + CQ) ] th  quantile 

of  F (D)  into  demand  D  under  risk  condition.  Let  q  =  C  / 

-•  _  o 

(C  +C  ),  then  we  have: 
s  o 

S*  =  F~ 1 (q)  . 

Although  we  don't  know  the  demand  distribution,  each  suc¬ 
cessive  decision  period  provides  another  observation  from  the 
demand  distribution,  and  we  will  assume  that  an  observation 
of  demand  is  indeed  obtained  each  period.  As  we  acquire  more 
data  we  might  be  able  to  estimate  S*  =  F  (q)  as  a  function 
of  the  given  data.  Two  kinds  of  approachs  are: 

1.  a  parameteric  approach  to  hypothesizing  the  form  of 
demand  distribution, 

2.  a  non-parametric  approach  to  find  the  quantile. 

For  the  first  case,  the  normal,  and  exponential  distribu¬ 
tions  have  been  found  to  be  of  considerable  value  in  des¬ 
cribing  demand  distributions.  The  normal  distribution  has 
been  found  to  describe  many  demand  functions  at  the  factory 
level;  the  exponential,  at  the  wholesale  and  retail  levels 
[3] .  Of  course,  these  distributions  should  not  be  auto¬ 
matically  applied  to  any  demand  distribution.  Statistical 
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tests  should  establish  the  basis  for  any  standard  assumption 
concerning  a  demand  function. 

We  are  directing  our  attention  to  the  case  where  the 
expected  value  solution  is  not  applicable  because  of  lack  of 
available  information.  Further,  if  we  had  enough  information 
for  a  fit  of  a  distribution,  we  should  have  a  good  estimate 
of  the  quantile  we  seek.  Hence,  the  second  case,  the  non- 
parametric  approach,  will  be  considered  to  find  the  optimal 
S*  and  switching  period. 
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III.  DEVELOPMENT  OF  THE  ESTIMATORS  ON 
THE  BASIS  OF  ORDER  STATISTICS 


As  mentioned  in  the  previous  chapter,  we  are  subject 
to  find  the  optimum  (minimum  cost  solution)  of  a  newsboy 
problem  without  knowing  the  q-th  quantile  of  the  demand  dis¬ 
tribution.  Consequently  a  non-parametric  or  distribution- 
free  approach  is  suggested  to  find  the  solution.  In  this 
chapter,  we  will  introduce  the  notion  of  order  statistics  as 
a  basis  for  estimating  this  quantile,  and  suggest  five  can¬ 
didate  quantile  estimators  which  may  be  used  in  newsboy 
solutions . 


A.  ESTIMATORS  USING  ORDER  STATISTICS 


The  non-parametric  approach  we  shall  follow  is  on  the 
basis  of  order  statistics  and  although  it  is  applicalbe  to 
both  continuous  and  discrete  random  variables,  we  shall  direct 
our  attention  to  the  continuous  case.  We  are  interested  in 
the  case  where  the  population  is  not  known,  but  we  have  a 
random  sample  of  size  n  from  its  unknown  cumulative  distribu¬ 
tion  function  F(X).  Let  X^,X2/...,X  be  a  random  sample. 

Then  X  ( 1)  <_  X(2)  <_  .  .  .  <_  X(n),  where  the  X(i)  are  arranged 
in  order  of  increasing  magnitudes,  and  are  defined  to  be  the 
order  statistics.  Let  X^  be  the  q-th  quantile  of  the  popu¬ 
lation.  The  estimate  of  X^  can  be  expressed  as  a  function  of 


x1,x2. 


.  ,x  , 

'  n 


i .  e . 
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Estimate  of  X 

q 


f  (X]_,X2,  . .  .  ,Xn)  . 


We  shall  consider  now  a  function  F  (x)  which  is  construc- 

n 


ted  from  the  sample  values  of  , • - . , XR.  The  function 

Fn(x)  is  called  empirical  cumulative  distribution  function. 

Let  x.,xn,...,x  denote  the  observed  values  of  X,,XV...,X  . 
l  2  n  1  z  n 

Then  for  each  number  x(-°°  <  x  <  °°)  ,  the  value  of  Fn(x)  is 

defined  to  be  the  proportion  of  observed  values  in  the  sample 

which  are  less  than  or  equal  to  x.  It  follows  from  the 

Glivenko-Cantelli  Theorem  that,  as  n  approaches  infinity  the 

empirical  cumulative  distribution  function  Fn(x)  converges 

uniformly  to  F(x).  Since  we  are  concerned  with  an  unknown 

cumulative  distribution  function  F(x)  which  is  continuous, 

some  type  of  smoothed  curve  of  Fn(x)  might  yield  a  reasonable 

estimator  of  F(x).  Let  q  be  the  [C  /(C  + C  ) ] th  quantile, 

o  s  o 

where  we  wish  to  find  the  estimate  of  X  so  that 


F(estimate  of  X^)  =  q.  This  is  represented  in  Figure  3. 

Now,  we  are  interested  in  the  estimator  which  gives  us  mini¬ 
mum  cost,  but  we  have  no  information  about  it  yet.  Let  us 
turn  our  attention  to  the  expected  cost  as  the  function  of  the 
estimate  of  X^.  Let  perturbation  be  the  difference  between 
the  estimate  of  X  and  true  value  of  X  .  Then  the  expected 

q  q 

cost  of  perturbation  is  proportional  to  the  variance  of  the 
given  unbiased  estimator  [4].  The  properties  of  unbiasedness 
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F  (x) 


Figure  3.  Estimate  of  X  on  the  Smoothed 
Curve  of  Empi?ical  CDF,  Fn(x) 

and  minimum  variance  should  be  considered  for  selecting 
estimators . 

B.  RATIONALE  FOR  THE  ESTIMATORS 

In  this  section,  we  shall  suggest  five  different  quantile 
estimators  on  the  basis  of  order  statistics.  In  a  statisti¬ 
cal  evaluation  of  these  estimators,  it  is  desirable  to  use 
the  uniform  (0,1)  distribution,  since  the  expectation  and 
variance  of  its  order  statistics  are  simple  and  exact. 
Investigating  these  properties  of  order  statistics  using  other 
distributions  is  not  easy..  If  we  have  a  distribution  similar 
to  the  uniform,  or  if  we  have  a  large  number  of  observations, 
some  approximations  may  be  possible.  Three  points  which  are 
near  each  other  can  be  assumed  to  be  linear  and  some  estimator 
using  order  statistics  would  be  approximately  unbiased. 
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We  shall  examine  the  properties  of  the  order  statistics 
from  the  uniform  (0,1)  distribution,  and  we  shall  consider 
three  neighboring  order  statistics.  Let  X(r)  '  an<^ 

X(r+i)  be  the  (r-l)th,  r-th,  and  (r+l)th  order  statistics 
respectively.  Define  the  value  r  as: 


r  =  [nq  +  0 . 5]  , 


(6) 


where  [X]  denotes  the  largest  integer  of  X.  Let  us  consider 
an  estimate  of  X  as  the  linear  combination  of  those; 

q 


(estimate  of  X  ) 

q 


aX 


(r-1) 


+  bX 


(r) 


+  cX 


(r+1) ' 


where  a,  b,  and  c  are  nonnegative  and  sum  to  unity.  There 
can  be  many  possible  ways  to  take  estimators,  but  we  shall 
examine  typical  ones  considering  the  unbiasedness  and  the 
magnitude  of  variance. 

A 

1.  The  First  Candidate  Estimator:  X 

-  q 

First  of  all,  the  estimator  should  be  unbiased, 

A 

i.e.,  E[X  ]  =  X  .  If  we  take  the  r-th  order  statistic, 

q  q 

X^  as  the  first  candidate  estimator,  i.e., 


(7) 


this  is  an  unbiased  estimator  since  E[X^^]  =  r/(n+l)  for 
the  r-th  order  statistic.  The  coefficients  are  a  =  c  =  0  and 
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b  =  1  in  this  case.  The  first  candidate  estimator  is  repre¬ 
sented  in  Figure  4. 


F  (x) 


2.  The  Second  and  Third  Candidate  Estimators:  X  ,  X 

- — - — - — -  q  q 

We  shall  propose  (r-l)th  and  (r+l)th  order  statistics 
as  the  second  and  third  candidate  estimators: 


Xq  "  X(r-1)  ' 


(8) 


and 


Xq  "  X  (r+1)  * 


(9) 


Even  though  these  are  not  unbiased  estimators,  we  are  not 
certain  which  one  gives  us  better  performance  since  we  are 
considering  the  properties  on  the  uniform  distribution  and 
there  is  some  round  up  error  for  taking  the  value  r. 
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3.  The  Fourth  Candidate  Estimator:  X 

-  q 

We  shall  propose  an  equally  weighted  linear  combina¬ 
tion  of  the  (r-l)th  and  (r+l)th  order  statistics  as  the 
fourth  candidate  estimator: 


X 


= 


(r-1) 


x(r+l) ] 


(10) 


This  estimator  is  an  unbiased  estimator,  since 


E[Xq] 


Ir£zi  +  £±1] 

2  ln+l  -  1 


n+1' 


r 

n+1  ' 


and  has  the  minimum  variance  using 


cov  [X 


(r) 


r  (n  -  s  -  1) 
(n+1) 2 (n+2) 


of  the  uniform  distribution  [5] .  The  variance  of  this  esti¬ 
mator  is  given  as  follows: 


Var [X  ] 

Si 


x .  ,  v  x,  , , 

Var  [' 1  "r~1  ~2+  (r+1)] 


(n-r) (2r-l) + (r-1) 
2 (n+1) 2 (n+2) 


This  estimator  is  represented  in  Figure  5. 

4.  The  Fifth  Candidate  Estimator:  X 

-  q 

We  shall  propose  an  equally  weighted  linear  combina¬ 
tion  of  the  (r-l)th,  r-th,  and  (r+l)th  order  statistics  as 
the  fifth  candidate  estimator,  or 


20 


F(x) 


3  [X(r-l) 


+  X 


(r) 


X(r+1)] 


(11) 


This  is  also  an  unbiased  estimator  and  the  variance  of  this 
estimator  is  given  by: 


Var  [X  ] 

q 


_  (n-r) (9r-4)  +  (5r-4) 

9 (n+1) 2 (n+2 ) 


The  difference  in  variance  between  the  fourth  estimator  and 

the  fifth  estimator  is  equal  to  ttt? — rrri — nvr*  When  n  =  1  and 

^  18 (n+1) (n+2) 

-4 

n  =  10  this  difference  is  approximately  0.09  and  4x10  , 

respectively.  The  difference  converges  rapidly  to  zero  as  n 
increases.  Since  the  quantity  difference  may  be  assumed  to 
be  negligible  and  the  variance  is  only  concerned  with  uni¬ 
form  distribution,  it  is  hard  to  say  which  estimator  gives 
us  better  performance. 
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In  summary,  although  there  may  be  many  possible 
candidate  estimators  depending  upon  different  forms  of 
linear  combination,  it  may  be  reasonable  to  examine  five 
previously  proposed  estimators.  They  can  be  shown  with 
properties  as: 


1) 

The 

First  Estimator: 

X 

q 

X^  (unbiased) 

2) 

The 

Second  Estimator: 

X' 

& 

II 

X '  ,  .  (biased) 

(r-1) 

3) 

The 

Third  Estimator: 

X*1 

& 

II 

X(r+1)  (biased) 

4) 

The 

Fourth  Estimator: 

X! 

II 

IIX(r-l)+X(r+l)> 

(unbiased  min.  variance) 

5) 

The 

Fifth  Estimator: 

X 

q 

TtX(r-l)+X(r)+X(r+l)) 

(unbiased) . 

Having  proposed  these  quantile  estimators,  it  is  of 
interest  to  see  how  well  they  perform  in  the  context  of 
newsboy  decision  making,  under  cost  criteria.  This  is  the 
subject  of  the  next  chapter. 
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IV.  EVALUATION  OF  THE  CANDIDATE  ESTIMATORS 


In  this  chapter,  the  five  candidate  estimators  will  be 
tested  and  investigated  for  their  performance  in  the  newsboy 
problem  for  several  demand  distributions.  As  performance 
measures,  two  kinds  of  costs,  namely  an  average  cumulative 
cost  and  a  mixed  cost  which  we  shall  introduce  later  on,  will 
be  obtained  during  the  simulation  and  then  the  performance 
of  the  candidate  estimator  will  be  evaluated  using  the  rela¬ 
tive  cost  deviation  to  the  ideal  cost. 

A.  DESIGN  FOR  SIMULATION 

A  newsboy  simulation  will  consist  of  50  replications, 
each  of  which  simulates  50  decision  periods  for  a  given 
demand  distribution.  The  quantile  value  (q  =  Co/(Cs  +  Cq)) 
for  each  distribution  is  varied  from  0.1  to  0.9  in  increments 
of  0.2. 

To  find  the  value  r  and  consequently  to  find  the  estimate 
of  X  ,  some  modification  of  the  method  of  determining  r  is 
necessary  during  the  first  and  the  last  few  periods,  since 
there  is  an  insufficient  amount  of  data  to  apply  the  candi¬ 
date  estimators.  For  example,  if  the  value  r  for  finding 
the  estimate  is  equal  to  zero,  then  X will  be  chosen  as 
the  estimate  of  X  ,  and  if  r  =  n,  then  the  estimate  of  X 

q  q 

will  be  X,  .  and  so  forth. 

(n) 

For  each  period,  the  average  cumulative  cost  and  the 
mixed  cost  are  computed  for  seven  different  cases  which  are 
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represented  as  follows: 


1)  Risk  case  (Ideal  solution) , 


2)  Uncertainty  case  using  the  minimax  rule, 


3)  A  case  using  the  estimator  X  , 


4)  A  case  using  the  estimator  X  , 


5)  A  case  using  the  estimator  X^, 


6)  A  case  using  the  estimator  X^,  and 


7)  A  case  using  the  estimator  X  . 

q 


The  cost  for  each  period  is  averaged  over  50  replications,  and 
then  the  average  cumulative  cost  at  period  k  is  computed  by 
cumulating  the  average  cost  from  the  first  period  to  k-th 
period.  The  mixed  cost  at  k-th  period  is  given  by  summing 
the  average  cost  using  the  minimax  rule  up  to  (k-l)th  period 
and  the  average  cost  using  an  estimator  from  the  k-th  period 
to  the  last  period. 


B.  DISTRIBUTION  AND  PARAMETERS 

As  mentioned  previously,  one  source  has  suggested  that 
the  normal  distribution  may  be  considered  as  the  demand  dis¬ 
tribution  at  the  factory  level  and  the  exponential  distribu¬ 
tion  may  be  considered  as  the  demand  distribution  at  the 
wholesale  and  retail  levels,  and  so  on.  Five  characterized 
forms  of  demand  distributions  including  symmetric,  right-skewed, 
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and  left-skewed  will  be  selected  to  test  the  performance  of 
the  candidate  estimators  for  generality,  and  are  shown  in 
Figure  6. 


f  (x) 


Figure  6.  Characterized  Forms  of  Demand  Distributions 

Demand  distributions  used  in  the  simulation  and  their 
parameters  are  shown  in  Table  I.  The  parameters  for  the 
normal  distribution  were  selected  so  that  negative  demand 
should  not  be  generated. 


Table  I.  Distributions  and  Their  Parameters 


Distribution 

Parameters 

Uniform 

The  Range  over  [0,30] 

Normal 

Mean  35  and  Variance  100 

Exponential 

Mean  20 

Gamma 

a  =  2  and  A  =  0 . 1 

Beta 

A  =  15  and  B  =  5 

25 


For  each  demand  distribution,  the  99.9th  percentile  of 

demand  except  uniform  distribution  will  be  adopted  as  D 

c  max 

to  apply  the  minimax  rule  under  uncertainty.  Quantile  values 
q  =  CQ/(Cs  +  C  )  will  be  0.1,  0.3,  0.5,  0.7,  and  0.9  for  each 
distribution. 

C.  ALGORITHM  FOR  SIMULATION 

The  program  steps  are  ordered  in  the  following  way: 

1.  Define  a  new  demand  distribution. 

2.  Assign  the  quantile  value  and  maximum  demand. 

3.  Compute  the  optimum  in  cases  of  uncertainty  and  risk. 

4.  Generate  three  random  deviates  from  the  specified 
demand  distribution,  representing  demand  for  the  first 
three  periods. 

5.  Obtain  the  order  statistics. 

6.  Compute  the  estimates  using  the  candidate  estimators. 

7.  Generate  a  random  deviate  for  the  specified  demand 
distribution  for  the  next  period,  compute  the  cost  for 
the  period  and  save  it. 

8.  If  the  period  is  beyond  the  final  period,  go  to  Step  9; 
otherwise  go  to  Step  5. 

9.  If  the  number  of  replications  is  greater  than  the  maxi¬ 
mum,  go  to  Step  10;  Otherwise  go  to  Step  4. 

10.  Find  the  average  cost  for  each  period  over  total 
replications . 

11.  Compute  the  average  cumulative  cost,  the  mixed  cost, 
and  their  relative  cost  deviation  to  the  risk  cost. 

Then  tabulate  the  results. 
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12.  If  quantile  value  is  greater  than  0.9,  go  to  Step  13; 
otherwise  go  to  Step  2. 

13.  If  we  want  to  try  another  demand  distribution  go  to 
Step  1;  otherwise  stop. 

D.  DISCUSSION  OF  THE  SIMULATION  RESULTS 

In  this  section,  the  simulation  results  will  be  displayed 
and  evaluated.  As  mentioned  previously,  the  simulation  run 
was  performed  over  50  replications,  each  of  which  consisted 
of  50  consecutive  demand  periods.  From  these  trials,  the 
average  cumulative  costs  and  the  mixed  costs  were  computed 
for  various  estimators. 

Table  II  shows  the  average  cumulative  cost  for  50  periods 
for  the  given  quantile  values  and  specified  demand  distribu¬ 
tion.  Errors  in  estimating  D  were  not  considered.  This 

max 

performance  of  quantile  estimation  may  be  compared  directly 
with  "best  case"  minimax  results.  Within  a  distribution  we 
may  conclude  one  estimator  is  better  than  the  other  but  it 
is  not  easy  to  conclude  for  all  distributions  using  Table  II. 

Now,  a  measure  of  effectiveness  (MOE)  to  evaluate  the 
performance  of  the  estimators  for  all  distributions  is  needed. 
Let  us  define  the  relative  cost  deviation  of  an  estimator  as 
follows : 

Relative  cost  deviation  _  (cost  using  estimator)  -  (risk  cost) 
(to  the  risk  cost)  (risk  cost) 
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Table  II 


Average  Cumulative  Costs  for  50 
Periods  and  50  Replications 


\r>esc 

Dist 1 n\ 

Estimator 

Quantile 

(q) 

0.1 

0.3 

0.5 

0.7 

0.9 

Risk 

5.498 

12.737 

15.038 

12.456 

5.310 

Uniform 

Minimax 

5.498 

12.737 

15.038 

12.456 

5.310 

(0,3) 

X 

5.778 

13.624 

16.468 

13.833 

6.290 

„q 

X 

5.846 

14.095 

17.397 

15.117 

7.702 

=q 

X 

5.782 

13.369 

15.935 

13.174 

7.253 

_q 

x„ 

5.672 

13.366 

16.195 

13.541 

7.011 

_q 

X 

q 

5.687 

13.376 

16.205 

13.580 

6.604 

Risk 

7.128 

14.015 

15.997 

13.805 

6.937 

Normal 

Minimax 

11.511 

19.943 

16.648 

15.399 

9.426 

(35,100) 

X 

8.272 

15.686 

17. 818 

15.120 

8.006 

_q 

X 

8.790 

16.870 

19.473 

16.892 

9.417 

~q 

X 

7.755 

14.724 

16.868 

14.577 

9.474 

_q 

X 

7.934 

15.288 

17.586 

15.010 

8.835 

_q 

X 

q 

7.944 

15.232 

17.502 

14.980 

8.329 

Risk 

7.444 

19.657 

27.541 

28.892 

18.642 

Exponen¬ 

Minimax 

11.705 

42.088 

60.824 

57.014 

26.512 

tial 
(Mean  = 

X 

_q 

7.958 

20.447 

28.589 

30.453 

20.906 

20) 

X 

7.949 

20.425 

28.853 

30.434 

21.114 

~q 

X 

a 

8.389 

21.623 

29.833 

33.840 

25.094 

X 

8.045 

20.576 

28.591 

30.746 

20.950 

=q 

X 

q 

8.046 

20.568 

28.547 

30.447 

20.072 

Risk 

6.714 

16.287 

21.141 

20.748 

12.496 

Gainma 

Minimax 

7.024 

23.094 

36.970 

37.240 

18.098 

( a  =  2 

X 

7.253 

16.963 

22.275 

22.141 

13.649 

X  =  0  •  X  ) 

_q 

X 

7.303 

17.180 

22.496 

22.442 

14.138 

=q 

15.932 

X 

7.641 

17.761 

22.861 

23.690 

_q 

X 

7.305 

17.067 

22.100 

21.992 

13.863 

-q 

X 

7.297 

17.043 

22.098 

21.964 

13.299 

q 
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Table  II  (Continued) 


^Desc 

Dist 

Estimator 

r~H 

• 

o 

0.3 

Quantile 

0.5 

(q) 

0.7 

0.9 

Risk 

0.071 

0.140 

0.154 

0.129 

0.059 

Beta 

Minimax 

0.260 

0.545 

0.520 

0.246 

0.059 

(A  =  15 

B  =  5) 

X 

.q 

0.079 

0.147 

0.162 

0.135 

0.064 

X 

=q 

0.079 

0.150 

0.167 

0.141 

0.070 

X 

_q 

0.082 

0.150 

0.163 

0.136 

0.082 

X 

=q 

0.078 

0.147 

0.161 

0.135 

0.071 

X 

q 

0.078 

0.147 

0.161 

0.135 

0.067 
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This  relative  cost  deviation  may  be  used  as  an  MOE  to  evalu¬ 
ate  the  performance  of  an  estimator.  The  best  performance 
was  given  by  the  risk  solution  during  the  simulation,  an 
expected  result  with  regard  to  theoretical  considerations. 

Table  III  exhibits  the  relative  cost  deviation  for 

quantile  values  from  0.1  to  0.9  for  each  demand  distribution 

investigated.  As  a  special  case,  if  we  pursue  the  uniform 

demand  distribution  in  Tables  II  and  III,  the  risk  and 

minimax  approach  would  match  each  other  because  the  risk 

cost  solution  X  such  that  F(X  )  =  C  /(C  +  CQ)  t*ie  same 

as  the  minimax  cost  solution  given  by  [C  /(C  +C  )]D 

^  2  o'  s  o  max 

whenever  the  maximum  demand  D  is  determined  exactly.  If 

max  2 

we  read  through  the  third  column  of  overall  distribution  in 
Table  III,  an  estimator,  X  shows  the  best  performance  among 
candidate  estimators  at  quantile  0.1.  Also,  X^  did  better 
at  quantile  0.9.  Excluding  these  two  points,  the  estimator 
Xg  gave  the  best  performance  at  all  other  quantiles.  Within 
the  uniform  distribution  the  minimax  approach  shows  the  best 
performance  (since  it  is  identical  to  the  risk  solution) , 
but  it  shows  the  worst  performance  in  the  other  distribution. 
Similarly,  we  can  say  which  estimator  represents  the  best 
performance  within  a  distribution  by  reading  through  row 
by  row,  and  which  estimator  shows  the  best  performance  at 
a  quantile  by  reading  through  column  by  column. 

We  are  interested  in  the  estimator  which  gives  us  the 
best  performance  for  all  quantiles  over  all  distributions. 
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Table  III 


Relative  Cost  Deviation  From  The  Average  Cumulative  Cost 


\Desc. 
Dist  'n\ 

Estimator 

0.1 

Quantile  i 
0.3  0.5 

(q) 

0.7 

0.9 

Overall 

Quantile 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Uniform 

Minimax 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

(0.30) 

X 

q 

0.0509 

0.0696 

0.0949 

0.1105 

0.1845 

0.5104 

^  u. 

X 

q 

0.0634 

0.1066 

0.1569 

0.2136 

0.4505 

0.9909 

X 

q 

0.0517 

0.0496 

0.0597 

0.0576 

0.3658 

0.5844 

_ M 

Xq 

0.0318 

0.0493 

0.0770 

0.0871 

0.3203 

0.5655 

X 

q 

0.0345 

0.0502 

0.0776 

0.0902 

0.2435 

0.4960 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Normal 

Miramax 

0.6149 

0.4229 

0.0405 

0.1154 

0.3589 

1.5526 

(35,100) 

X 

q 

0.1605 

0.1192 

0.1136 

0.1962 

0.1542 

0.6427 

X 

~q 

0.2331 

0.2037 

0.2173 

0.2236 

0.3576 

1.2353 

X 

a 

0.0879 

0.0505 

0.0544 

0.0559 

0.3658 

0.6145 

X 

q 

0.1130 

0.0908 

0.0994 

0.0782 

0.2736 

0.6640 

— M 

X 

q 

0.1145 

0.0868 

0.0941 

0.0850 

0.2008 

0.5812 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Exponen¬ 

Minimax 

0.5724 

1.1411 

1.2085 

0.9733 

0.4222 

4.3175 

tial 
(Mean  = 

xq 

0.0690 

0.0402 

0.0381 

0.0540 

0.1215 

0.3228 

20) 

X 

a 

0.0678 

0.0391 

0.0476 

0.0603 

0.1326 

0.3474 

X 

Q 

0.1269 

0.1000 

0.0832 

0.1713 

0.3461 

0.8275 

X 

rr 

0.0806 

0.0467 

0.0381 

0.0642 

0.1238 

0.3534 

=4 

X 

q 

0.0809 

0.0463 

0.0365 

0.0538 

0.0767 

0.2942 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Gamma 

Minimax 

0.0412 

0.4180 

0.7488 

0.7949 

0.4483 

2.4562 

(a  =  2 

X  =  0.1) 

X 

q 

0.0803 

0.0415 

0.0536 

0.0671 

0.0922 

0.3347 

X 

rr 

0.0878 

0.0548 

0.0641 

0.0817 

0.1314 

0.4198 

X 

rr 

0.1381 

0.0905 

0.0814 

0.1418 

0.2749 

0.7267 

__4 

X 

n 

0.0881 

0.0479 

0.0454 

0.10600 

0.1094 

0.3508 

—4 

X 

q 

0.0869 

0.0464 

0.0453 

0.0586 

0.0642 

0.3014 
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Table  III  (Continued) 


\Desc. 

Dist'n'\ 

Estimator 

0.1 

Quantile 
0.3  0.5 

(q) 

0.7 

0.9 

Overall 

Quantile 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Beta 

Minimax 

2.6322 

2.9002 

2.3714 

0.9108 

0.0085 

8.8291 

(A=  15 

B=  5) 

X 

0.1023 

0.0519 

0.0549 

0.0481 

0.0864 

0.3436 

X 

0.1123 

0.0753 

0.0877 

0.0952 

0.1823 

0.5528 

X 

0.1481 

0.0736 

0.0558 

0.0575 

0.3852 

0.7202 

_q 

X 

_q 

0.0916 

0.0508 

0.0459 

0.0468 

0.2019 

0.4368 

X 

q 

0.0917 

0.0481 

0.0462 

0.0451 

0.1307 

0.3618 

Risk 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Overall 

Minimax 

3.8657 

4.8822 

4.3752 

2.7944 

1.2379 

17.1554 

Dist 'ns 

X 

q 

0.4630 

3.3224 

0.3551 

0.3749 

0.6388 

2.1542 

^  tl 

X 

a 

0.5644 

0.4795 

0.5736 

0.6744 

1.2543 

3.5462 

X 

q 

0.5527 

0.3642 

0.3345 

0.4841 

1.7378 

3.4733 

_ tl 

X 

q 

0.4049 

0.2855 

0.3058 

0.3453 

1.0290 

2.3705 

— -u. 

X 

q 

0.4085 

0.2771 

0.2997 

0.3327 

0.7159 

2.0346 
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One  way  to  compare  is  to  note  that  the  performance  of  an 
estimator  overall  quantiles  over  all  distributions  can  be 
represented  as  the  total  sum  of  the  relative  cost  deviations. 
Looking  at  the  last  column  for  overall  quantiles  over  all 
distributions  it  may  be  concluded  that  the  best  performance 
(minimum  sum)  is  given  by  the  estimator  X^  =  +X(r) 

+  X(r+]jl  for  aH  quantiles  over  all  distributions  as  shown 
in  Table  III,  because  the  minimum  relative  cost  deviation  is 
the  closest  to  the  risk  cost  (ideal  cost  solution) .  For 
specific  quantiles,  of  course,  occasionally  other  estimators 
performed  better.  In  an  application  of  the  newsboy  problem 
under  initial  uncertainty,  one  has  information  about  the 
quantile  q  but  often  not  about  the  form  of  the  distribution. 
Accordingly,  the  information  at  the  foot  of  Table  III  may 
provide  information  in  selecting  an  estimator. 

E.  DISCUSSION  OF  SWITCHING  RULES 

When  solving  a  newsboy  problem  repeatedly  under  uncer¬ 
tainty  conditions  we  should  probably  use  the  minimax  rule  for 
the  first  period  since  there  is  no  applicable  amount  of  data 
to  use  in  quantile  estimation.  Each  successive  decision 
period  provides  an  additional  piece  of  demand  data  which 
will  contribute  to  quantile  estimation.  Therefore  we  are 
interested  in  the  number  of  decision  periods  that  should 

elapse  before  we  apply  the  proposed  estimator  X  in  place  of 

Si 

the  minimax  rule.  From  the  simulation  results.  Table  IV 
shows  the  mixed  cost  with  switching  at  period  25,  the  middle 
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Table  IV 


The  Mixed  Cost  for  the  25-th  Period 


\Desc. 

DistT"\ 

Estimator 

0.1 

0.3 

Quantile 

0.5 

(q) 

0.7 

0.9 

Risk 

274.879 

636.863 

751.893 

622.810 

265.514 

Uniform 

Minimax 

274.879 

636.863 

751.893 

622.810 

265.514 

(0.30) 

280.729 

643.679 

761.937 

634.742 

272.186 

.V 

X 

_q 

282.951 

648.889 

766.665 

640.032 

277.903 

X 

q 

280.185 

641.266 

759.937 

631.729 

269.095 

X 

q 

279.242 

643.131 

760.970 

633.978 

271.334 

— >1 

X 

q 

279.595 

643.109 

761.187 

634.082 

271.463 

Risk 

356.404 

700.766 

799.844 

690.275 

346.833 

Normal 

Minimax 

575.550 

997.144 

832.229 

769.934 

471.301 

(35,100) 

A 

X 

q 

461.583 

828.894 

820.863 

738.497 

413.906 

X 

q 

477.707 

835.620 

825.044 

742.866 

420.868 

Z  ^ 

X 

q 

455.163 

825.489 

819.134 

736.412 

410.975 

i 

X 

q 

460.505 

828.510 

820.079 

737.839 

412.314 

460.535 

828.517 

820.258 

737.909 

412.595 

Risk 

372.210 

982.875 

1377.038 

1444.599 

932.085 

Exponen¬ 

Minimax 

585.271 

2104.407 

3041.196 

2850.687 

1325.620 

tial 

X 

458.956 

1453.137 

2073.317 

2028.197 

1118.620 

(Mean  = 

20) 

^q 

X 

CT 

458.897 

1453.506 

2075.067 

2027.887 

1115.982 

~  M 

X 

a 

461.103 

1456.723 

2076.325 

2036.363 

1146.552 

X 

a 

458.524 

1453.449 

2072.446 

2027.401 

1121.038 

— 

X 

q 

458.556 

1453.220 

2072.413 

2027.286 

1119.067 

Risk 

335.684 

814.349 

1057.038 

1037.387 

624.824 

Gamma 

Minimax 

351.179 

1154.714 

1848.499 

1861.989 

904.908 

(a  =  2 

X 

344.370 

983.554 

1414.183 

1400.872 

759.497 

A  =  0.1) 

^q 

X 

Q 

344.799 

985.234 

1415.761 

1402.143 

759.026 

X 

Q 

348.497 

983.682 

1410.905 

1406.584 

768.769 

X 

344.654 

982.802 

1410.740 

1400.742 

754.786 

=q 

Xq 

344.504 

982.942 

1411.690 

1400.479 

755.480 
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Table  IV  (Continued) 


\Desc. 

Distts\ 

Estimator 

0.1 

0.3 

Quantile 

0.5 

(q) 

0.7 

0.9 

Risk 

3.573 

6.993 

7.699 

6.437 

2.948 

Beta 

Minimax 

12.988 

27.275 

26.003 

12.299 

2.973 

(A=  15 

B=  5) 

X 

q 

7.691 

15.796 

15.779 

9.190 

3.021 

X 

=q 

7.726 

15.789 

15.779 

9.250 

3.056 

X 

_q 

7.723 

15.840 

15.809 

9.214 

3.015 

J3 

7.680 

15.788 

15.781 

9.191 

3.014 

X 

q 

7.678 

15.789 

15.779 

9.190 

3.015 
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period  in  our  fifty-period  simulation.  The  mixed  costs 
here  are  the  sum  of  costs  for  the  first  24  periods  given  by 
the  minimax  rule  and  the  costs  given  by  an  estimator  from 
the  25-th  to  the  last  period.  When  we  examine  various  dis¬ 
tributions  in  Table  IV  the  minimum  mixed  cost  is  obtained 
by  the  risk  approach,  which  is  to  be  expected,  and  the 
proposed  estimator,  X  ,  appears  to  provide  the  minimum 

q 

mixed  cost  among  candidate  estimators,  although  by  a  narrow 
margin. 

Now,  we  wish  to  find  the  number  of  periods  to  retain 

the  minimax  rule  before  applying  the  proposed  estimator  X  . 

Si 

Mixed  costs  were  computed  for  various  switching  periods  to 
find  the  period  which  gives  the  minimum  mixed  cost  for  each 
quantile  of  every  distribution  except  uniform  distribution, 
since  the  proposed  estimator  X^  doesn't  give  better  results 
than  the  minimax  rule  for  the  uniform  distribution.  Table 
V  shows  the  best  switching  periods  (minimax  to  estimator  X  ) 
for  the  quantiles  and  distributions  used  in  the  simulation 
during  40  periods. 


Table  V 

Best  Switching  Periods  Using  the  Proposed  Estimator  X 


Distribution 

0.1 

0.3 

Quantile 

0.5 

(q) 

0.7 

0.9 

Normal 

4 

7 

27 

11 

11 

Exponential 

6 

4 

4 

4 

4 

Gamma 

13 

5 

4 

4 

6 

Beta 

4 

4 

4 

4 

48 
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Looking  at  Table  V,  we  see  that  the  best  switching  period 
is  27  at  quantile  0.5  of  the  normal  distribution.  This 
result  is  caused  by  the  fact  that  the  minimax  rule  behaves 
well  around  the  median  when  demand  has  a  symmetric  distri¬ 
bution.  The  occurrence  of  switching  period  48  at  quantile 
0.9  of  the  beta  distribution  can  be  explained  by  the  fact 
that  with  the  beta  distribution  used  in  the  simulation, 
the  optimal  S*  obtained  by  the  minimax  rule  is  close  to 
the  ideal  S*  at  this  quantile,  i.e.,  the  minimax  S*  at  0.9 
is  equal  to  0.8767  and  the  ideal  S*  is  0.8660. 

We  shall  suggest  6  periods,  4,  5,  6,  7,  11,  and  13  as 
the  candidates  for  the  switching  period  for  all  quantiles 
over  all  distributions.  Later  periods  are  not  included  in 
the  candidates  since  the  estimator  X  shows  much  better 

q 

performance  (with  reference  to  the  foot  of  Table  III  for 
nearly  all  quantiles  over  all  distributions) .  Table  IV 
shows  the  mixed  costs  obtained  by  the  estimator  for  the 

q 

25th  period.  Similar  tables  (not  shown)  were  computed  for 
various  other  periods,  and  showed  that  the  mixed  costs 
increased  as  the  period  increased.  We  propose  the  relative 
mixed  cost  deviation  for  a  candidate  period  using  the  esti¬ 
mator  X  as  an  MOE .  It  may  be  reasonable  to  conclude  that 

q 

the  period  giving  the  minimum  is  the  most  preferable.  Table 
VI  shows  the  relative  mixed  cost  deviation  at  each  candidate 
period  for  all  quantiles  over  all  distributions  using  the 
estimator  X  . 

q 
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Table  VI 


Relative  Mixed  Cost  Deviation  Using  the  Estimator 


Switching 

Period 

Distribution 

Estimator 

0.1 

0.3 

Quantile  (q) 
0.5  0.7 

0.9 

Total 

Uniform 

.0345 

.0502 

.0776 

.0902 

.2435 

.4960 

Normal 

.1145 

.0868 

.0941 

.0880 

.2008 

.5812 

4 

Exponential 

X 

q 

.0809 

.0463 

.0365 

.0538 

.0767 

.2942 

Gamma 

.0869 

.0464 

.0453 

.0586 

.0642 

.3014 

Beta 

.0917 

.0481 

.0462 

.0451 

1307 

.3618 

Overall 

.4085 

.2778 

.3327 

.3327 

.7159 

2.0346 

Uniform 

.0324 

.0471 

.0675 

.0762 

.1720 

.3952 

Normal 

.1201 

.0868 

.0751 

.0673 

.1462 

.4955 

5 

Exponential 

X 

q 

.0708 

.0548 

.0561 

.0669 

.0851 

.3337 

Gamma 

.0576 

.0422 

.0555 

.0713 

.0708 

.2971 

Beta 

.1235 

.1003 

.0902 

.0570 

.0868 

.4578 

Overall 

.4044 

.3312 

.3444 

.3387 

.5606 

1.9793 

Uniform 

.0317 

.0457 

.0626 

.0693 

.1229 

.3322 

Normal 

.1255 

.0863 

.0663 

.0644 

.1145 

.4570 

6 

Exponential 

X 

q 

.0657 

.0697 

.0742 

.0765 

.0894 

.3755 

iTVitytta 

0477 

0774 

31  nl 

•  Uu  /  J 

•  V  Oi't 

Beta 

.1715 

.1591 

.1400 

.0790 

.0608 

.6104 

Overall 

.4435 

.4035 

.4106 

.3706 

.4650 

2.0932 

Uniform 

.0300 

.0420 

.0547 

.0601 

.1020 

.2888 

Normal 

.1311 

.0846 

.0533 

.0572 

.1084 

.4346 

7 

Exponential 

X 

q 

.0738 

.0924 

.0980 

.0944 

.0987 

.4573  . 

Gamma 

.0423 

.0487 

.0754 

.0924 

.0847 

.3435 

Beta 

.2306 

.2237 

.1935 

.0973 

.0481 

.7932 

Overall 

.5078 

.4914 

.4749 

.4014 

.4419 

2.3174 
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Table  VI  (Continued) 


Switching 

Period 

Distribution 

Estimator 

0.1 

0.3 

Quantile  (q) 
0.5  0.7 

0 . 9  Total 

Uniform 

.0262 

.0323 

.0384 

.0422 

.0502  .1891 

Normal 

.1727 

.0952 

.0368 

.0526 

.0988  .4561 

11 

Exponential 

X 

q 

.0859 

.1641 

.1779 

.1594 

.1292  .7165 

Gairma 

.0194 

.0797 

.1253 

.1384 

.1041  .4669 

Beta 

.4119 

.4456 

.3793 

.1705 

.0375  1.4448 

Overall 

.7161 

.8169 

.7577 

.5631 

.4196  3.2734 

Uniform 

.0247 

.0280 

.0323 

.0383 

.0446  .1679 

Normal 

.1916 

.1047 

.0322 

.0544 

.1070  .4899 

13 

Exponential 

X 

q 

.1057 

.2066 

.2228 

.1932 

.1413  .8696 

Gairma 

.0910 

.0930 

.1502 

.1627 

.1158  .6127 

Beta 

.5236 

.5666 

.4823 

.2166 

.0362  1.8253 

Overall 

.9366 

.9989 

.9198 

.6652 

.4449  3.9654 
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The  relative  mixed  cost  deviation  at  each  candidate 

period  is  shown  in  Table  VII.  The  minimum  deviation  is 

given  at  period  five  and  the  deviation  is  increasing  as  the 

period  increases  from  the  fifth  period.  Using  this  result 

we  may  propose  that  the  switching  period  from  the  minimax 

rule  to  applying  the  estimator  X  should  be  the  fifth  period. 

Si 


Table  VII 

The  Relative  Mixed  Cost  Deviation 
At  Each  Candidate  Period 


Period 

Cost  Deviation 

4 

2.0346 

5 

1.9793  (minimum) 

6 

2.0932 

7 

2.3174 

11 

3.2734 

13 

3.9654 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 


The  minimax  approach  is  a  well  known  approach  to  decision 
making  under  uncertainty.  For  a  newsboy- type  problem,  we 
developed  five  candidate  estimators  which  could  be  used  as 
data  aggregates,  and  evaluated  them  by  comparing  with  the 
minimax  rule.  These  candidate  estimators  were  tested  for 
quantiles,  0.1,  0.3,  0.5,  0.7,  and  0.9  of  several  distribu¬ 
tions:  uniform,  normal  (symmetric),  exponential,  gamma 

(right-skewed)  and  beta  (left-skewed) .  Since  all  possible 
quantiles  and  other  kinds  of  distributions  were  not  covered, 
there  may  be  some  amount  of  loss  of  generality. 

With  these  limitations,  as  discussed  in  Chapter  IV,  the 

estimator  X  =  i[X,  ..  +X,  .  +X,  , .  .  ]  developed  on  the  basis 

q  3  (r-1 )  (r)  (r+1)  v 

of  order  statistics  showed  the  best  performance  among  those 
considered.  The  simulation  results  suggest  that  if  a  single 
decision  procedure  over  all  quantiles  is  needed,  the  esti¬ 
mator  X  is  recommended  and  the  time  for  switching  to  this 
estimator  should  be  at  the  fifth  period.  The  proposed 
procedure  may  be  summarized  in  the  following  way: 

(1)  Estimate  D  . 

in  3.x 

(2)  Compute  the  quantile  q  =  Co/(Cs  +  Cq). 

(3)  Apply  the  minimax  rule  S*  =  <2*Dmax  for  the  first 
four  periods. 

(4)  Compute  the  value  r  such  that  r  =  [nq  +  0.5]  from  the 
fifth  period  where  n  is  the  period  at  which  we  wish 
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to  estimate  the  demand  and  [X]  denotes  the  largest 


integer  of  X. 


(5)  Compute  the  order  statistics  X(r)'  an<^  x 

and  apply  the  proposed  estimator  X  : 

q 


(r+1) 


X 


q 


+  x 


(r) 


+  X 


(r+1) 


]  • 


In  this  study,  it  was  implicitly  assumed  that  the  value 
of  demand  for  each  period  was  obtainable  even  if  D  >  S.  For 
further  study,  the  case  when  demand  data  greater  than  the 
quantity  S  is  not  available  is  worthy  of  investigation. 

Also,  a  non-parametric  approach  which  can  handle  extreme 
quantiles  (e.g.,  0.01  or  0.99)  and  extremely  skewed  distri¬ 
butions  could  be  investigated.  The  comparison  between  a 
parametric  approach  when  the  demand  data  are  assumed  to  fit 
a  probability  distribution  and  non-parametric  approach  is 
also  recommended.  Finally,  future  simulation  work  on  this 
problem  should  also  record  cost  variances,  so  that  a  statis¬ 
tical  analysis  could  be  made  of  the  simulation  results. 
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APPENDIX 


PROGRAM  LISTING  OF  THE  SIMULATION 
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